Critical conductance of two-dimensional chiral systems with random magnetic flux 
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The zero temperature transport properties of two-dimensional lattice systems with static random 
magnetic flux per plaquette and zero mean are investigated numerically. We study the localiza- 
tion properties and the two-terminal conductance and its dependence on energy, sample size, and 
magnetic flux strength. The influence of boundary conditions and of the oddness of the number 
of sites in the transverse direction is also studied. For very long strips of flnite width, we flnd a 
diverging localization length in the middle of the energy band at _E = and determine its critical 
exponent v = 0.35±0.03. A previously proposed crossover from a power-law to a logarithmic energy 
dependence can be excluded from our data, at least for energies |_E| > 10~^°. For square systems, 
the sample averaged scale independent critical conductance {qc) turns out to be a function of the 
amplitude of the flux fluctuations whereas the variance of the respective conductance distributions 
appears to be universal. We flnd a critical conductance (gc) — 1.49e^//i for the strongest possible 
disorder. 

PACS numbers: 73.23.-b, 71.30.+h, 72.10.-d 
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I. INTRODUCTION 

The transport properties of charged quantum particles 
in two-dimensional systems with various types of disorder 
are of considerable interest in a variety of experimental 
and theoretical situations. In particular, the presence of 
a static random magnetic flux with zero mean has been 
of much concern recently in connection with bond disor- 
dered Anderson models with either real or complex hop- 
ping terms fii^i^i^i^i^iii^i^ii^ with the composite- fermion 
picture of the fractional quantum Hall effect at half- 
filling,— the critical behavior at the quantum phase 
transition of spin-split Landau levelsfi^ and with the 
gauge field theory of high- Tc superconductivityJ^ii^ In 
addition, far-reaching relations between the low energy 
chiral limit of a quantum chromodynamic (QCD) parti- 
tion function and a large N limit of the random matrix 
theorjii^iil as well as between the electrical conductance 
in disordered media and spontaneous chiral symmetry 
breaking in QCD have recently become apparent J^ii^ 

Concerning the random flux model, there exists an 
extensive list of valuable contributions to this intricate 
problem (see, e.g., Ref. |2(| and references therein), but 
a definite picture started to emerge only recently, at 
least for quasi-ld (QID) samples,^*' Results for true two- 
dimensional systems are scarce and precise numerical es- 
timates are still missing. A consensus has been reached 
on the notion that all electronic states are localized for 
such systems where in addition to the random magnetic 
flux also random diagonal disorder is present j^i^i.^^ How- 
ever, in the absence of diagonal disorder, it has also been 
shown that the random flux model with Gaussian dis- 
tributed and (5-correlated magnetic fields can be mapped 
onto a nonlinear a model of unitary symmetry so that all 
electronic states should be localized)^ The recognition 
of a special chiral symmetry that can emerge in systems 
with an underlying bi-partite lattice, so that the eigen- 
values appear in pairs ±6^)^ has considerably augmented 



our view of the possible situations a random flux model 
can assume,2^i22iS^iSa 

Our aim is to investigate a lattice model with static 
random magnetic fluxes and to numerically calculate the 
two-terminal conductance and the localization properties 
for energies close to the band center. We want to study 
the role of the chiral symmetry and to clarify the possible 
dependence on boundary conditions (BC). In addition, 
we address the influence of an odd or even number of 
lattice sites. For QID systems, we will check the asser- 
tion that the Lyapunov exponents do not come in pairs 
for samples with an even widthj^ We will also look for 
the crossover proposed for the energy dependence of the 
localization lengthi^ Finally, we calculate the size de- 
pendence of the conductance of square systems and show 
that at the band center the conductance converges to the 
critical value {gc) — lA9e^/h. 

The paper is organized as follows. In Section IIIII we 
study the spectrum of Lyapunov exponents (LE) in the 
quasi-ld limit and for square samples and discuss how the 
physical symmetry of the system depends on the bound- 
ary conditions, the parity of the width of the system, and 
energy of the electron. In Section [TVl we find that E — 
is a critical point: the smallest LE does not depend on 
the width L of the lattice. For L odd and Dirichlet BC, 
we also calculate the critical exponent for the divergence 
of the localization length of the two-dimensional system. 
We find a value ly — 0.35 ± 0.03 which is close to the crit- 
ical exponent for the Anderson bond disordered model'^'— 
and also in agreement with the one obtained with a dif- 
ferent method for the random-flux model which has been 
reported recently.^ For QID systems of finite width L, 
the localization length diverges as £_ cx i|ln(|£;|Li/'')|. 
In Section |V] we present our data for the critical two- 
terminal conductance. Although the scale independent 
mean value (gc) depends on the strength of the magnetic 
field fiuctuations /, the variance of the corresponding 
distributions Pc{9,f) turns out to be universal. We also 
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confirm the unusual length dependence of the mean con- 
ductance for systems with L odd and DBC.^° Concluding 
remarks are given in Section IVII 

II. MODEL AND METHOD 

The two-dimensional (2d) motion of non-interacting 
particles subject to a perpendicular random magnetic 
field is described by a Hamiltonian 

^ = ~ {^x{c\^+a^C.jy^ + Cm-g^C^) (1) 

m 

with nearest neighbor hopping, defined on the sites m of 
a 2d square lattice, where the width L (a;-direction) and 
the length (z-direction) of the sample is measured in 
units of the lattice constant a. The value of the hopping 
term in the x-direction is i^: = 1 if not stated otherwise, 
and the energy is given in units of = 1. The oper- 
ators and create or annihilate a Fermi particle 
at site TO, respectively. The complex hopping terms are 
chosen such that the magnetic flux (in units of the flux 
quantum h/e) through an individual plaquette is given 
by the sum of the random Peierls phases along the two 
bonds in the z-direction 2tt4>„i = am,m+a, - oim+a^ ,m+a, ■ 
The random fluxes are distributed uniformly according to 
-//2 < 0m < //2, where < / < 1, with probabihty 
density p{<j>m) = 1// so that its second moment is /^/12, 
and the average magnetic flux through the system is zero. 
The randomness is maximal for f — I. 

Without attached leads, the model H]) exhibits chiral 
unitary symmetry for Dirichlet boundary conditions in 
both directions. The chirality is destroyed when periodic 
boundary condition are imposed along any direction pro- 
vided the number of sites in this direction is odd. TableU 
summarizes the various situations. The chiral symmetry 
is always broken by an additional on-site disorder. 

In the following, we study numerically the quantum 
transport of electrons with energy E through the 2d sys- 
tem defined by the Hamiltonian ([T]). For a given length 
of the system L^, we calculate the dimensionless two- 
terminal conductance via the relational 

^ 1 

g-Tr{rtr} = ^ (2) 
^ cosh (Xi/2) 

In Eq. ([2|), T is the transmission matrix and the Xi pa- 
rameterize its eigenvalues. The electrons propagate in 
the z-direction and N is the number of open channels. 
Dirichlet (DBC) or periodic (PBC) boundary conditions 
are imposed in the transversal direction. In the limit of 
Lz/L oo, the parameters Xi converge to the quantities 
Zi X {Lz/ L)^^^ where Zi is the zth Lyapunov exponent 
(LE) which characterizes the exponential decrease of the 
wave function of quasi-ld systems. Oseledec proved'^'* the 
convergence for the eigenvalues of the transfer matrix. 
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TABLE I: The symmetries of the model Hamiltonian ([T]) de- 
pend on the boundary conditions and on the oddness of the 
number of sites. In the absence of leads, the possible sym- 
metry classes are unitary (U), chiral unitary (CU), and chi- 
ral unitary with an extra eigenvalue that appears at i5 = 

(r \ I I N 20.27.28.29 



Zi = Yimi^^/i^^^ Zi{Lz/ L). For sufficiently large L^/L, 
the Zi{Lz/L) are self-averaging quantities. The smallest 
positive LE zi is related to the localization length and 
represents the key parameter of finite-size scalingj^i^^ 

Since the calculation of the transmission probability re- 
quires two semi-infinite (ideal in our case) leads attached 
to the left and right of the sample, the boundary con- 
dition in the propagation (z) direction are neither PBC 
nor DBC. We expect, however, that the boundary condi- 
tions in the transversal (a;) direction affect the transport 
properties of the system considerably. 

Our data, both for the conductance and for the pa- 
rameters Xi, support the conjecture that (i) the system 
possesses chiral unitary symmetry only at the band cen- 
ter £;= for DBC, and for PBC with L even. The 
chirality of the E = state is confirmed by our data for 
the parameters Xi. In particular, we find that the prob- 
abihty p{x{) does not decrease to zero when xi — > 0. We 
will later discuss that this behavior is typical for the chi- 
ral symmetry class, {ii) There exists a critical point at 
the band center for L odd and DBC. Since this critical 
point is due to the chiral symmetry of the model, we ex- 
pect the criticality also for L even. This expectation is 
supported by our numerical data for the smallest LE zi . 
For L odd and PBC, the critical state at E = Q should 
disappear due to the unitary symmetry. 



III. LYAPUNOV EXPONENTS 

Figure [T] shows the spectrum of Lyapunov exponents 
\zi\ for quasi-ld systems with Dirichlet and periodic BC 
in the transverse direction and with either odd {L — 65) 
or even (L = 64) system width. For L even, the spectrum 
is degenerate at the band center for both Dirichlet and 
periodic BC 

\z2^-l\ = \z2^\ = cx[i-^/2] (Lcvcn). (3) 
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FIG. 1: The energy dependence of the spectrum of Lyapunov 
exponents \zi\ of the transfer matrix. Dirichlet and periodic 
BC are imposed in the transversal direction, and f ~ 1- Left: 
L — 65, right: L = 64. Dashed lines indicates the values of 
the Lyapunov exponents for E = 0. Note that 2:1 = for L 
odd and Dirichlet BC. 



For L odd, we obtain at the band center that 
( c X Int [i/2] {L odd, DBC) 

\^^\^{ (4) 

[ c/2 X [i - 1/2] {L odd, PBC). 

From Fig. [T] and later from Fig. [8] we see that c « 2.68. 
As is shown in Fig. [1] the degeneracy is removed for non- 
zero energy. In the transfer matrix method, we calculate 
only positive Lyapunov exponents. Since the LE appear 
in pairs, we have also doubly degenerate LE {—Zi,—Zi) 
in the negative part of the spectra. 

While the form of the spectrum of LE for odd L and 
PBC is typical for unitary symmetry;^ the degeneracy 
of the spectra, observed in all three other cases, indicates 
chiral symmetry i^^i^'^i^^ 

For DBC, the chirality is confirmed also by the anal- 
ysis of the distribution of parameters Xi, calculated for 
finite length L^- Since we are able to calculate only the 
absolute value of the LE, we cannot distinguish from the 
present data whether the system possesses chiral unitary 
(CU) or unitary (U) symmetry. Fortunately, we can es- 
timate the physical symmetry from the analysis of the 
distribution of the parameter Xi, calculated for systems 
of finite length L^. 

As discussed in Refs. [13 and [H, for weak disorder 
the probability distribution p{{x}) is determined by the 
Dorokhov-Mello-Pereyra-Kumar equation^^i^S 
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FIG. 2: The probability distribution p(|2:i|) for E = and 
L = 33, 65, and 129 (data points). Solid lines are p{xi) 
for L = 65 and S = 10"*, 10"^, and 10~^ (from the left). 
The last distribution is compared with the Wigner surmise 
{xi)Wi{xi) = -Isexp — ^s^, where s = x\/{xi). The in- 
set shows p[xi) for E = 10"^ and L = 33, 65, and 129. 
Shown is also the distribution p{x2 — xi) for E = 10""^ and 
L = 129 which is almost identical with the Wigner surmise 
{a;i)W2(a;i) = -^-s^exp— for unitary ensemble. 



Here, £ is the mean free path and J is the Jacobian 



J = 



Uky , |sinh(a;j - XkW 



n 



(CU) 



I sinh 



sinh^^fe|2n,|sin(2a;,)| (U). 



(6) 

The main consequence of the absence of the repulsion 
term sin(2a;i) in the Jacobian ^ is that the spectrum 
of Xi spans over the entire real axis: the Xi can be both 
positive and negative when the system possesses chiral 
unitary symmetry. In the ordinary unitary systems, all 
values of positive, being reflected from the ori- 

gin by an additional term in the Jacobian. Clearly, in 
the case of unitary symmetry, p{xi) — > when xi — > 0, 
but p{xi — 0) is non-zero in the case of chiral symme- 
try. Since we are not able to calculate the sign of the 
parameters xi for a given sample, we plot in Fig. [2] the 
distribution of the absolute value \xi\. Dirichlet BC are 
imposed in the transversal direction. For E ^ the 
distribution does not depend on the system size. If xi 
possesses both positive and negative values, the distribu- 
tion p{xi) is Gaussian with a mean value = 0. This 
agrees with our data for the quasi-ld systems where we 
find zi = 0. Therefore, we conclude that the system 
possesses chiral symmetry. 

However, the form of the distribution p{xi) changes 
qualitatively when the energy differs from zero. As is 
shown in Fig. [21 already for E — 10"^ the distribution 
p(|a;i|) decreases to zero when |a;i| — > 0. This confirms 
that the Jacobian given by Eq. ([6]) contains also the repul- 
sion term cx sin(2a::). Consequently, the system changes 
the symmetry from chiral unitary to unitary and all pa- 
rameters Xi become positive. As shown in Fig. [21 the 
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FIG. 3: The probability distribution p{xi) for a system with 
PBC at J5 = 0. The system size is 65 x 66 and U = 0.99. 
p{xi) agrees well with the Wigner surmise for orthogonal and 
p{xi — X2) (shown in the inset) for unitary ensembles.— 
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FIG. 4: The probability distribution pi2, defined by Eq. ((Tj, 
for is = (shaded area) and for E — 0.001. The size of the 
system is 66 x 66. Dirichlet BC are used in the transversal 
direction. Shown are also distributions p{xi) and p{x2) for 
E = 0.001. Note that p ^ when x ^ 0. This confirms that 
the system possesses different physical symmetry for E = 
and i5 7^ 0, in agreement witb^S. 



distribution p{xi) converges to the Wigner surmise Wi 
when either E or L increases. Also, the distribution of 
differences X2 — xi converges to the Wigner surmise W2- 
This behavior of ^(xi) and p{x2 — xi) is typical for the 
unitary universality classi^ 

In case of L odd and PBC in the transverse direction, 
the symmetry changes to unitary and the critical point 
at E = disappears. The transfer matrix algorithm does 
not enable us to calculate the parameters Xi for E = 
and PBC due to the kz = eigenmode of the transfer ma- 
trix in unperturbed leads. This mode disappears either 
when E ^ or when an anisotropy in the hopping terms 
is applied. Using a small anisotropy in the a;-direction, 
tx = 0.99, we confirmed that the statistics of p(2;i) and 
p{x2 — xi) follow the Wigner surmises also at the band 
center. Fig. [3] shows the respective distributions p{xi) 
to be Wi and p{x2 — xi) is W2- This is in contrast to 
the situation with DBC where the distribution changes 
qualitatively on approaching E — 0. 

In case of L even, the analysis is more difficult since 
we expect the mean values of the first two parameters 
xi and X2 to have the same absolute value but with an 
opposite sign. So, we cannot distinguish between |xi| and 
\x2\ in our analysis of a given sample. To overcome this 
problem, we calculate for iVstat realizations the common 
probability distribution. 



Pl2ix) = —— J2 -^(^ - l^ll) + - l^2|), (7) 
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therefore that the distribution 
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NstsA 



P12 



(^) = ^E'^(^-^i) + '^(-^-^2) (8) 



is of the form 



^12(2:) 



e 2<T -(_ e 2<T 



(9) 



This expectation is confirmed also by Fig.[5l which shows 




of the parameters |a;i| and \x2\ for square system L x 
L with L — 66 and DBC in the transversal direction 
(Fig. [5]). We see that the probability pi2{x) (shown by 
the shaded area) is non-zero when x ^ 0. We expect 



FIG. 5: The distribution pi2{x) for systems with even width 
L = 66. The length of the system is Lz = 66 (top), Lz = 132 
(middle) and Lz = 1000 (bottom). The solid line is the fit of 
Pi2(x), given by Eq. 
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E, a^lEIL 

FIG. 6: The energy dependence of the spectrum of Lyapunov 
exponent zi of the transfer matrix. Hard wall boundary con- 
ditions are imposed in the transversal direction and / = 1. 
The width of the system is given in the legend. The data 
scale to the universal curve on the right hand side, described 
by Eq. (I13|l . with critical exponent v — 0.342 and 02 — 0.05. 
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FIG. 7: The energy dependence of the spectrum of Lyapunov 
exponents zi of the transfer matrix. Hard wall boundary 
conditions are imposed in the transversal direction and the 
flux strength is / = 0.5. The width of the system is given in 
the legend. The data scale to the universal curve (rhs) given 
by Eq. (HSJ with u = 0.359 and 02 = 0.30. 



how the probability distribution changes when the sys- 
tem length increases. 

For square systems, we obtain the distribution shown 
in Fig.[4]which, in the limit of L^/L — > 00 transforms into 
two Gaussian peaks. For longer systems, p{xi,X2) devel- 
ops into two isolated Gaussian peaks centered around the 
mean values, (a;i) = — (X2). In analogy to the odd L case, 
a non-zero energy breaks the chiral symmetry also in the 
even L situation. A similar statistics was observed also 
for PBC with small anisotropy (not shown), which con- 
firms the existence of the chiral symmetry also for L even 
and PBC. 

We conclude that the random flux model with Dirichlet 
EC possesses at the band center E = a chiral unitary 
symmetry. The spectrum of the Lyapunov exponents is 
given by the relations 

Zi=cx (-1)'+^ Int [{i + l)/2] (L even) (10) 

and 

z, = cx {-lyint [i/2] (i odd, DBG), (11) 
in agreement with previous theoretical considerations.-^' 

IV. CRITICAL REGIME AND EXPONENT 

Since zi = for i? = 0, i odd and DEC in the 
transversal direction, the system is in the critical regime 
with a diverging correlation length 



To estimate the critical exponent v, we calculate Zi as a 
function of energy E and of the system width L. We ex- 
pect, in agreement with the single parameter scaling j'^^i^° 
that zi is a function of the ratio L/£^{E) only. As is shown 
in Fig. [HI all numerical data can be fitted by the universal 
function 



zi(E,L) 



ai 



\\n{a2\E\Ly-)y 



(13) 



with three fitting parameters Oi, 02, and v. From the 
scaling analysis we observed that 



0.35 ±0.03. 



(14) 



^^\E\ 



(2d). 



(12) 



More detailed information of the analysis is presented 
in Table [TTl To estimate the accuracy of our result, we 
repeated the scaling analysis with reduced input data 
sets. A similar value of the critical exponent was obtained 
also for systems with weaker magnetic flux disorder / — 
0.5 (see Fig. [T]). Due to the smaller values of zi, we have 
to simulate much longer quasi-ld systems in order to get 
data with reasonable accuracy. 

Although the calculated values of for / — 0.5 dif- 
fer slightly from those obtained for / = 1.0, we do not 
interprete this difference as a non-universality of the crit- 
ical exponent.^ Rather we assume that this difference is 
due to the limited accuracy of our numerical data and/or 
fitting procedure. Indeed, as shown in Table |TT1 the es- 
timated value of the critical exponent depends on the 
choice of the input ensemble defined by zi^ax and £min, 
and decreases slightly when larger values of zi are ex- 
cluded. 

We did not find any crossover from a power-law to 
a more complicated E dependence of the localization 
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33 
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4/36 


33 


0.20 


0.384 


11/42 


33 


0.25 
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26/45 



TABLE II: Numerical estimate of the critical exponent v for 
two different strengths of the random flux amplitudes, / = 1.0 
and / = 0.5. Only data for L > Lmin and with z\ < zimax 
are considered in the scaling analysis. -Fmin is obtained from 
the minimum of the fitting function, A^'data is the number of 
data. The accuracy of the critical exponent in each fitting 
procedure is of the order of 10~^. 
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FIG. 8: The smallest Lyapunov exponent z\ as a function of 
the system width L for L even (top) and L odd (bottom). The 
data confirm that z\ does not depend on the system width. 
This either implicates the existence of a critical point at the 
band center in all three cases or a finite size effect due to the 
limited system size. 



length for 2d as proposed in Ref. [13 and discussed in 
Refs. H and H. Our numerical data cannot be fitted to 
the one parameter scaling function zi{E,L) = zi(L/£) 
with a localization length ^ cx exp ^J\n{EQ/E)^ Since 
we analyze a very narrow energy interval around the band 
center (as small as \E\ ~ 10~^°), we do not expect that 
the crossover from the observed power-law to the pro- 
posed logarithmic energy dependence of the localization 
length exists in our situation. 

Since zi determines the localization length of the quasi- 
Id system, ^ = 2/zi , we see from Eq. (|13|) that the local- 
ization length diverges as 

U{E)o^Lx\\n{a2\E\L^/'')\ (QID), £; ^ (15) 

for a given system width L. A logarithmic divergence is 
typical for Anderson bond disordered models.— 

While the existence of the critical state at = for L 
even is commonly accepted)^ we do not expect the same 
for L odd and PBC, since the system possesses unitary 
symmetry in this case. To describe the property of the 
E = Q state, we plot in Fig. [8] the L dependence of the 
smallest LE zi for L even (both Dirichlet and periodic 
BC) and L odd (PBC). In all three situations, we do not 
observe any L dependence of the smallest LE zi. We 
believe that this indicates that the localization lengths 
considerably exceeds the available system sizes so that 
no final conclusions can be reached. 

The scaling analysis is very difficult in the case of L 
even. As is shown in Fig. [HI the L-dependence of zi 
is highly non-trivial for non-zero energies in accordance 
with Ref. 3. The scaling seems to work only in the limit 
of L ^ 00. We disagree on the observation^ that the LEs 
at i? = do not come in degenerate pairs. In contrast 



we find the difference between the two LEs to be smaller 
than the accuracy of our calculations. 

The criticality of the E = state for the Dirichlet BC 
will be supported also by the size dependence of the mean 
conductance, discussed in the next Section. 

V. CONDUCTANCE 

Figure [TOl shows the size dependence of the sample av- 
eraged critical conductance {gc) for square systems Lx L 
{L odd) and three values of the randomness strength 
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FIG. 9: The L dependence of the smallest Lyapunov exponent 
zi for L even and for various values of the energy E. The non- 
monotonous L dependence disables the scaling analysis in this 
case. Our data are consistent with Fig. 1 of Ref. [3. 
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FIG. 10; The critical value of the mean conductance (qc) 
as a function of the system's size L x L with L odd (open 
symbols) and L even (full symbols) at iJ = with DBC, and 
various strengths of the random field /. The data show that 
the critical conductance does not depend on the parity of L, 
but depends on /. For completeness, we add also data for 
E = 0.0001 (L odd) and E = 0.001 (L even) to show that the 
conductance decreases with L when E ^ 0, indicating that 
the system is in the localized regime in the limit of L cx). 



/. The energy is £^ = and Dirichlet BC are consid- 
ered. Our data confirm that (g) converges to an L in- 
dependent critical value (gc) which, however, does de- 
pend on the strength of the randomness /. For the 
largest possible disorder / = 1 we obtain {gc) = 1.49, 
a value larger than the 2d symplectic case.^'^ Figure [TT] 
shows the / dependence of the mean conductance for 
squares of size 257 x 257. It also shows that the variance 



var g 



(g) is universal and independent on /. 
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FIG. 11: The critical value of the mean conductance ((?c) for 
square samples 257 x 257 at £ = and various values of the 
random field /. Dirichlet BC are used in the transversal direc- 
tion. The inset shows the probability distribution p{g — (g)) 
for three different values of /. The width of the distribution 
varg ^ 0.187 does not depend on /. 
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FIG. 12: The length dependence of the mean conductance (g) 
for QID systems. The systems width is L = 65, the energy 
E = 0. One clearly sees the crossover from the 1/L (dotted 
line) to l/v^Lr(dashed line) dependence predicted by Ref . [20l. 



We find a value var 17 « 0.187 which is in agreement with 
those obtained earlier by Ohtsuki et al^ and Furusakii^ 
We observe, however, an increase of var g for very small 
/, which can be explained by finite size effects due to a 
large mean free path. 

We also plot in Fig.[TO]the size dependence of the mean 
conductance for squares with even L. Within the ob- 
tained accuracy, {gc) does neither depend on the parity 
of L nor on the boundary conditions in agreement with 
Ref. iS). Contrary to the band center, the conductance de- 
creases always with increasing system size whenever the 
energy lies outside the band center. 

We also analyzed the length dependence of the mean 
conductance {g) and of the mean of the logarithm of the 
conductance (In g) for systems with hard wall transversal 
boundary conditions and E — Q. Since {xi) = in this 
case, we expect that the system possesses an infinite lo- 
calization length also in the quasi-ld limit;^ Therefore, 
the mean conductance {g) should not decrease exponen- 
tially when the system length increases. 

Our results shown in Fig. [12] confirm the relations pre- 
dicted theoretically^'^ 



{g) = MIL, 



and 



{g) = ^iM/i-KL,) 
in the limit of {g) « 1 and {g) <C 1, respectively. 



(16) 



(17) 



VI. SUMMARY 

We investigated two-dimensional electron systems with 
static random magnetic flux and showed numerically that 



8 



the transport properties depend on the parity of the sys- 
tem's width L and on the transverse boundary condi- 
tions. For Dirichlet boundary conditions, we confirmed 
by the analysis of the statistical properties of the quan- 
tities a;, which parameterize the eigenvalues of the trans- 
mission matrix, that the system possesses chiral unitary 
symmetry at the band center. The chirality exists in case 
of Dirichlet boundary conditions for both L odd and even, 
and for periodic BC for L even only. But the chirality 
is always broken when the energy of the electron is non- 
zero. 

In case of chiral unitary symmetry, the 2d system with 
random magnetic flux possesses a critical point at the 
band center. We found that the localization length di- 
verges cx \E\~''^ when £^ ^ and calculated the criti- 
cal exponent v w 0.35 for L odd and Dirichlet BC. Our 
data do not confirm the existence of the crossover from 
the power-law to the logarithmic energy dependence of ^ 
predicted by Reflsol 



We also calculated the critical conductance of 2d sys- 
tems. At the band center, the mean conductance {g) con- 
verges to a size-independent critical value for both L odd 
and L even. Although the critical conductance does de- 
pend on the strength of the randomness, the fluctuations 
of the conductance appear to be universal. For non-zero 
energy, the mean conductance decreases with the system 
size, indicating a localized regime. Finally, for the quasi- 
Id systems with odd system width and Dirichlet BC, we 
confirmed the non-trivial length dependence of the mean 
conductance, proposed theoretically in Ref. [l^- 
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